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Abstract 

Cosmic ray electrons and positrons constitute an important component of the background for imaging atmospheric Cherenkov 
Telescope Systems with very low energy thresholds. As the primary energy of electrons and positrons decreases, their contribution 
to the background trigger rate dominates over protons, at least in terms of differential rates against actual energies. After event 
reconstruction, this contribution might become comparable to the proton background at energies of the order of few GeV. It is well 
Q_|known that the flux of low energy charged particles is suppressed by the Earth's magnetic field. This effect strongly depends on 
the geographical location, the direction of incidence of the charged particle and its mass. Therefore, the geomagnetic field can 
contribute to diminish the rate of the electrons and positrons detected by a given array of Cherenkov Telescopes. 

In this work we study the propagation of low energy primary electrons in the Earth's magnetic field by using the backtracking 
I— —I technique. We use a more realistic geomagnetic field model than the one used in previous calculations. We consider some sites 
relevant for new generations of imaging atmospheric Cherenkov Telescopes. We also study in detail the case of 5@5, a proposed 
low energy Cherenkov Telescope array. 
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1. Introduction 

\ The great success achieved by current systems of Imaging 
Atmospheric Cherenkov Telescopes (lACTs) has led ground- 
based gamma-ray astronomy to enter a period of great develop- 
ment, with the next generation of instruments being designed 
|to reach unprecedented sensitivities and angular resolution in 
a more extended energy range, particularly with lower energy 
thresholds {E,h). One of the main motivations for the construc- 
tion of lACTs with low is the study of distant extragalactic 
sources. Gamma-rays from such sources are attenuated by their 
■interactions with the background radiation present in the inter- 
galactic medium. Due to this effect, the spectrum is severely 
suppressed, e.g. the energy cutoff for an extragalactic object 
with redshift z - I can be as low as ~ 50 GeV llj]. For a gen- 
eral discussion of the physical motivations to lower the energy 
threshold of lACTs see, e.g. |2']. 

Current systems of lACTs like HESS, MAGIC and VER- 
ITAS have energy thresholds of ~ 100 GeV with very large 
collection area. For a customized trigger system, the energy 
threshold of MAGIC can be as low as ~ 25 GeV. At these ener- 
gies gamma-rays can also be detected by instruments on satel- 
lites. The Fermi LAT gamma-ray detector is able to detect pho- 
tons in the energy range from ~ 20 Me V up to more than ~ 300 
GeV Ist]. Its large field of view makes it a very eflicient detec- 
tor for the discovery of new sources. However, the collection 
area in satellites is very limited so the sensitivity worsens very 



rapidly above tens of GeV. Thus, although the energy range of 
Fermi overlaps with ground-based lACTs, the sensitivity at en- 
ergies of few GeV is low for both detectors. 

The most advanced project considering a low £,/, is the 
Cherenkov Telescope Array (CTA), the largest international ef- 
fort for the next generation of lACTs with an order of magni- 
tude better sensitivity than current systems 101 ■ In particular, E,h 
for CTA was planned to be ~ 10 GeV, although a more reahstic 
value for the present design is ~ 20 GeV. There are also specific 
proposals whose main technical objective is to lower the en- 
ergy threshold, hke STEREO ARRAY U, ECO- 1000 H and 
5@5 101 ■ For the 5@5 array it is shown that an £,/, in the range 
3-5 GeV can be achieved by 5 lACTs placed at 5 km of alti- 
tude. To test the feasibility of measurements at high altitudes, 
the OMEGA project is also under consideration ||8|. 

Cosmic ray protons and electrons are the most important 
background sources for the discrimination of gamma-ray show- 
ers developed in the Earth's atmosphere. Electron initiated 
showers are practically indistinguishable from those initiated 
by gamma-rays and thus their importance for gamma-ray as- 
tronomy. For energies below ~ 20 GeV, it is predicted that elec- 
trons could be differentiated from gamma-ray initiated showers 
101 . Protons are almost two orders of magnitude more numer- 
ous than cosmic electrons at energies below ~ 10 GeV, increas- 
ing toward higher energies as the electron spectrum is steeper 
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than the one corresponding to protons IlIOllllll . However, pro- 
tons need quite a bit more energy to produce the same amount 
of Cherenkov light, specially at the lowest energies, going from 
a factor ~ 5 (at ~ 1 TeV) to -60 (at few GeV) Sill. Addition- 
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ally, the trigger system of lACTs discriminate against protons, 
although this is less effective at the lowest energies due to fluc- 
tuations in the shower development and a drop in image inten- 
sity. In the end, for the lowest energies, electrons become the 
dominant component of the differential trigger rate for a low- 
energy-threshold lACT system ^ Q [lH. All proton events 
passing the analysis cuts are considered as gamma-rays, and 
thus their reconstructed energy will be lower than the true en- 
ergy. This causes the electron dominance to diminish, making 
the reconstructed flux of both electrons and protons comparable 
at the lowest energies ijst]. 

Cosmic ray electrons and positron^ constitute ~ 1 % of the 
total cosmic ray flux arriving at Earth in the GeV-TeV energy 
range. It is believed that the high energy component of the elec- 
tron flux is directly produced by galactic sources such as su- 
pernova remnants and pulsars lll4ll . Electrons can also be pro- 
duced by interactions of cosmic ray protons or light nuclei with 
the interstellar medium gas. Electrons undergo energy losses 
during their propagation in the interstellar medium. The main 
processes are; synchrotron radiation in the galactic magnetic 
field, inverse Compton scattering with photons from stars and 
the cosmic microwave background, bremsstrahlung with inter- 
stellar matter, and ionization. For energies grater than ~ 10 
GeV, the electron flux is dominated by the local environment 
because the attenuation length is reduced to the kpc scale lll411 . 

Depending on the location on Earth, the geomagnetic field 
can act as a shield for charged particles, suppressing the flux of 
low energy electrons and, in this way, diminishing their con- 
tribution to the background for lACTs. A first study of the 
geomagnetic field effects on the cosmic electron rate was per- 
formed by Cortina and Gonzalez lITsll . in the context of the 
MAGIC telescopes, for several geographical positions around 
the world by using the dipolar approximation of the magnetic 
field of the Earth. A better description of the geomagnetic field 
is provided by the International Geomagnetic Reference Field 

(iGRF) m, which is given as a multipole expansion up to lOth 

or 13th order, depending on the version under consideration. In 
this case, the problem cannot be solved analytically, and thus 
numerical methods are used. In particular, the backtracking 
technique is used to find the allowed and forbidden trajectories 

The shielding of charged particles is not the only effect 
caused by the magnetic field of the Earth at a given location. 
Extensive air showers develop in the atmosphere generating 
negatively and positively charged particles, particularly elec- 
trons and positrons being the most important for lACTs. These 
particles are deflected in opposite directions by the component 
of the geomagnetic field normal to the shower axis. This spread 
causes an additional dispersion in Cherenkov images recorded 
by lACTs and, consequently, diminishing the background sep- 
aration efficiency. There are extensive studies in the literature 
about this effect (e.g. ifisll ). which depends not only on the 
location on Earth but also on the telescope pointing direction. 



While this and other effects might be of more importance than 
the shielding effect, they are not considered in this paper 

In this work we study the suppression of the cosmic electron 
flux due to geomagnetic field effects for locations in the South- 
em hemisphere. Moderate to high altitude sites with clear skies 
are available only in South Africa and South America. Thus, we 
consider here three candidate sites in the southern hemisphere 
(which are being considered for the installation of CTA): El 
Leoncito (31:47 S, 69:28 W), San Antonio de los Cobres (SAC, 
(23:50 S, 66:16 W)), both in Argentina UJl, and the HESS site 
in Namibia ll20ll . We also study in detail the case of 5@5, for 
which we consider the sites Llano de Chajnantor in the Ata- 
cama desert, northern Chile |i7|, and SAC which is very close 
to the latter (less than 200 km southwest). 

2. The cosmic electron flux 

The electron flux has been measured by several experiments 
(see j 10;] for a compilation of experimental data). The most re- 
cent measurements with better statistics are from Fermi LAT 
I in and PAMELA (l^. Fermi LAT covers the energy range 
from 7 GeV to 1 TeV, whereas the corresponding PAMELA 
range is from 1 GeV to 0.625 TeV. Although above ~ 10 GeV 
the PAMELA spectrum seems to be softer than Fermi LAT's, 
they are consistent within uncertainties. For energies below 10 
GeV the data from older experiments fall within the range of 
the PAMELA and HEAT results. The latter can be taken as 
the extreme minimum case for the measured electron flux arriv- 
ing at Earth. For energies below ~ 10 GeV, the discrepancies on 
the electron flux measured by different experiments can be ex- 
plained, in part, by solar modulation effects suffered by incident 
electrons. 

In this work, the flux measurements from Fermi LAT, 
PAMELA and the low energy part (£ < 10 GeV) of HEAT 
are considered. Figure [T] shows the experimental data, where 
the error bars indicate statistical plus systematic uncertainties. 
For the numerical calculation an analytical expression of the 
electron and positron flux (/ - Jeie + Jpos) is used, which is ob- 
tained by fitting the data shown in figure [1] with the following 
function. 
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Hereafter, electrons will refer to both electrons and positrons. Any partic- 
ular case will be specified explicitly. 



where E is in GeV, J is in m^^s^^ sr^^GeV^^ and both 
{a,b,c,d] and {0o,y) sets of parameters are not independent, 
but related by the conditions that make the flux and its deriva- 
tive continuous at £ = 7 GeV. 

Figure [T] shows both fits considered here as the maximum 
and minimum cosmic electron fluxes: the first corresponding to 
the combination of PAMELA and Fermi LAT data (Jpfl- solid 
line), and the second to the combination of Fermi LAT and the 
low energy part of HEAT data {Jhfl- dashed line). Table [T| 
summarizes all parameters resulting from the fits, as specified in 
equation ([T). It is considered that any other measured flux falls 
in between these two fits so that all possibilities are covered 
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Figure 1: Electron flux multiplied by obtained by Fermi LAT, PAMELA and the low energy part (E < 10 GeV) of HEAT. Error bars include statistical and 
systematic uncertainties. The sohd line corresponds to the fit of data from PAMELA and Fermi LAT. The dashed line coixesponds to the fit of data from HEAT and 
Fermi LAT. 
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JpFiiE) 


1.12439 


-22.2423 


142.483 
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256.183 


3.1213 


Jhfl{E) 


0.42764 


-12.1735 
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-113.49 


208.317 
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Table 1: Parameters resulting from the fits of equation (U] for PAMELA & Fermi LAT (Jpfl) and for HEAT & Fermi LAT (Jhfl)- 



within uncertainties. The relevance of the goodness of these 
fits is discussed in next section. Nevertheless, to evaluate the 
geomagnetic field effects at a given site, the fit Jpfi is finally 
chosen as it is the more recent and better measured electron flux 
available at present. 

Although at energies of the order of a few GeV the flux is 
dominated by the electron component, the positrons can con- 
tribute in a non negligible way to the lACTs background. Fig- 
ure |2] shows the positron fraction, 6{E) = Jpos{E)/{Jp„s{E) + 
JeieiE)) as a function of primary energy, obtained by PAMELA 
121 and Fermi LAT Q. Also shown is a fit of the PAMELA 



data with the function, 

log 6(E) = po + /5i log £■ + p2 log^ E, 



(2) 



where E is in GeV, po = -1.078, pi = -0.542 and p2 = 0.352. 
Note that the fit is consistent with the Fermi LAT data and, 
therefore, it is reliable up to 200 GeV, the maximum energy 
reached by Fermi LAT. 

3. Numerical technique and results 

Charged particles propagating through the interstellar 
medium towards the Earth are deflected by the geomagnetic 
field. The deflections suffered by these particles depend on their 
rigidity, which is defined as: R - pcjZe, where c is the speed of 
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Figure 2: Positron fraction as a function of primary energy measured by 
PAMELA and Fermi LAT. The solid line corresponds to a fit of the data with 
the function in Eq. (2)- 



light, p is the momentum and Ze is the charge of the particle. In 
this way, particles coming from a given direction of incidence 
with rigidities smaller than a given value are prevented to reach 
the Earth's surface. This minimum value of rigidity is called 
rigidity cutoff for that particular location and direction. 

The dominant component of the geomagnetic field is origi- 
nated in the interior of the Earth. As a first approximation it can 
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be modeled by a dipole located in the center of the Earth. In 
this case, it is poss ible to find an analytical expression for the 
rigidity cutoff !25\. A better approximation is to consider an 
eccentric dipole for which it is also possible to solve the prob- 
lem analytically. However, precise calculations require the use 
of more sophisticated models, like the one provided by the In- 
ternational Geomagnetic Reference Field (IGRF) lfl6ll . In this 
case, the problem cannot be solved analytically and then nu- 
merical methods are used to obtain the trajectories of particles 
UtIi . In this work, the propagation of cosmic electrons in the 
magnetic field of the Earth is performed by using the program 
rjI95 lE6ll written by Smart and Shea. The trajectories are cal- 
culated by using the backtracking technique, i.e. the trajectory 
of an electron, arriving at a given geographical location, with 
energy Eq, and arrival direction «, is calculated by propagating 
a positron with initial energy £0 and direction -n. Similarly, 
the trajectories of positrons are calculated by propagating elec- 
trons. The geomagnetic field implemented in the program cor- 
responds to a multipole expansion up to IQth order (the external 
magnetic field, e.g. solar, is not included) and the equations of 
motion are solved by using the fifth order Runge-Kutta method. 
The electrons and positrons are propagated from an altitude of 
20 km. 

As an example of allowed and forbidden rigidities (or en- 
ergies) for a particular location and direction, figure |3] shows 
the trajectories corresponding to electrons falling vertically at 
El Leoncito. The red line corresponds to electrons with energy 
E^ie - 9 GeV, the magenta to 10.65 GeV, and the blue to 12 
GeV. As can be seen from the figure, all energies but the one 
corresponding to Eeie = 9 GeV are allowed. 




Figure 3: Trajectories corresponding to electrons falling vertically at El 
Leoncito with energies 9 GeV (red line), 10.65 GeV (magenta line), and 12 
GeV (blue line). 

For an electron arriving to a given location on Earth with pri- 
mary energy E, zenith angle 6 and azimuth angle (p, the func- 
tion Teie{E, 9, <p) is defined to take the value 1 if the trajectory 
is allowed, and if it is forbidden. Similarly, the function 
TposiE, 6, (p) is defined for positrons. 

Figure|4]shows T^e as a function of energy for electrons, also 
falling vertically at El Leoncito. The figure clearly shows the 
penumbra region, which is defined as the energy transition be- 
tween the forbidden and allowed energy regions. In the exam- 



ple, the penumbra is centered at about 1 1 GeV and has a width 
of w 0.8 GeV, delimited by two energy values. 



Ef and Ef, 



such that r^/e is zero below E^^ and one above E"^*" 



o 

II 

-e- 

0.6 

O 

II 

<Xl 

td 0.4 



9.S 



10.5 11 11.5 

E [GeV] 



Figure 4: 
Leoncito. 



Tgig as a function of energy for vertical electrons falling at El 



The function T, and also Ei and Eh, strongly depends on 
the geographical location of the impact point and the direction 
of incidence. Figure |5] shows contour plots in AitofF projection 
of E"^" as a function of azimuth and zenith angles, for all sites 
under consideration. The maximum zenith angle considered 
here is 60°. The azimuth angle is measured clockwise looking 
down from the North. The figure shows the well known East- 
West eff'ect, i.e. for negatively charged particles arriving with 
the same zenith angle, the cutoff (or in this case Eff) is smaller 
for the East direction. For El Leoncito E^'' ranges from ~ 8 
GeV to ~ 25 GeV, for SAC from -8.5 GeV to ~ 26.8 GeV and 
for Namibia from ~ 5.2 GeV to ~ 11 GeV. Therefore, the largest 
(i.e. the best) values of Eh correspond to SAC. However, the 
difference with El Leoncito is very small. On average, £'^'' for 
Namibia is about a factor two (or even more) smaller than the 
value corresponding to the other two sites. 

Similar results are obtained for positrons. In this case the ar- 
rival directions with the smallest values of E'^" are in the West 



region, as can be seen from figure |ATT]( Appendix A for SAC. 

For a Cherenkov telescope with low enough energy thresh- 
old, the expected detection rate of electrons for a given site 
strongly depends on the arrival direction when geomagnetic 
field effects are included. In order to quantify this effect, the 
following parameter is introduced, 

F(E,h, e, 4>) = FMh, 6, 4,) + Fpos(E,h, 0, 4>), (3) 

where, 

dET,i,{E,e,cl))J,u{E) 



Fele{Eth,0, 



F pos(Eih, 0, 



dE J(E) 



J^~'dETpUE,0,<p)Jpos(E) 
C""" dE J(E) 

•JLti, 



(4) 



(5) 



F{E,h, 9, (p) is the ratio between the detection rate of electrons 
with and without including geomagnetic field effects, for an 
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Figure 5: Contour plots of E'^ as a function of azimuth and zenith angles for 
the three sites under consideration. The azimuth angle is measured clockwise 
looking down from the North. 



ideal detector with detection area A{K) - Ao &{E - E,k), 
where Aq is a constant and 0(x) is the Heaviside function (i.e. 
0(x) = 1 if X > 0, and 0(x) = if < 0). Here is 
the maximum energy considered. Note that the fluxes for elec- 
trons and positrons can be written in terms of both the total flux 
and the positron fraction, i.e. JeidE) = (1 - 6{E)) J{E) and 
Jpos(E) = 6(E) J(E). 

In order to calculate F, equation ([1) is used together with 
the functions T^ie and Tpos calculated with the rjI95 program. 
The maximum energy used to compute F is E,„„y- - 200 GeV, 
the upper limit of the energy interval for the positron fraction 
reported by Fermi LAT (see Fig. |2}- 

To show an example the lowest proposed energy threshold is 
taken, £,/, = 3 GeV Figure |6] shows contour plots in Aitoff' 
projection of F, as a function of azimuth and zenith angles, for 
the electron flux Jpfi (see tablelTJ. As expected, smaller values 



of F are obtained for El Leoncito and SAC. This is due to the 
larger values of obtained for these two sites, in comparison 
with Namibia. Although the difference between F for the first 
two sites is small, on average, F is smaller for SAC. Similar 
results are obtained by using the fitted electron flux Jhfl (see 
tablelB. 
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Figure 6: Contour plots of f as a function of azimuth and zenith angles for 
E,i, = 3 GeV. Black crosses indicate the directions of incidence for the maxi- 
mum and minimum values of F, obtained for 6,„ax = 40° and for 6,„nv = 60° . 
The electron flux used in the calculation con'esponds to the one obtained by 
fitting the PAMELA and Fermi LAT data ( Jpfl)- 



For the characterization of a given location a study of F as a 
function of £,/, will be presented below. For that study, the min- 
imum and maximum values of F are taken as representative of 
each location, this is F„,i„ and F,„a v. These quantities will cor- 
respond to two fixed positions on the sky (i.e. to given azimuth 
and zenith angles) which depend on the maximum zenith an- 
gle under consideration. The lACT technique is regularly used 
for observations at zenith angles smaller than ~45°. For larger 
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zenith angles there is a gain in collection area and a deteriora- 
tion of the image quality. It is not obvious how these two effects 
change the telescope performance so simulations are needed 
1 271 12811 . In any case, it is usual for lACTs to observe with 
zenith angles below ~60°. Therefore, two cases are considered 
in this work in order to obtain both the maximum and minimum 
values of F for each site: < 40° and 6 < 60°. These points 
are indicated by crosses in the cases shown in figure |6] for both 
values of d„,ax- 

For £,/, - 3 GeV the reduction on the arriving electron rate 
can be significant. In table |2] the minimum and maximum frac- 
tions of arriving electrons as defined in equation (l3) are summa- 
rized for all locations under study, using both electron fluxes, 
JpfL and Jhfl, and for 0„,„i = 40°. As expected, all values 
found using Jhfl are larger than those obtained with Jpfi as 
the latter is the largest electron flux considered, with the differ- 
ences in table |2]being in the range 15 - 24%. With these result, 
it is clear that to improve the fit goodness for Jpfi and Jhfl 
is not relevant. Even if a much better fit is done, particularly 
for Jhfl in the region of £ ~ 7 GeV, the effect on F would 
not be significant (estimated as <5%). Similar conclusions are 
reached when O^a.^ = 60° is considered. 
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F 

^ mm 


^ max 




JpFL 


Jhfl 


JpFL 


Jhfl 


El Leoncito 


0.031 


0.037 


0.12 


0.14 


SAC 


0.026 


0.031 


0.10 


0.12 


Namibia 


0.14 


0.16 


0.33 


0.38 



Table 2: Fraction of aniving electrons, F, with energies above 3 GeV (i.e. £,;, = 
3 GeV). Minimum and maximum values of F are considered for 6„,nv = 40° 
(see text) and for both electron fluxes, Jpfl and Jhfl- 




Figure 7: Contour plots of FposlF^ig as a function of azimuth and zenith angles 
for SAC, and using the electron flux Jpfl- 



these two cases. Figure [8] shows F as a function of energy 
threshold for the electron flux Jpfl- Solid lines correspond to 
Omax - 40° and dashed lines to 6„,„x - 60°. From the same 
figure it can be seen that, if the energy threshold of a given 
Cherenkov telescope is larger than ~ 27 GeV and d,„ax = 60°, 
F - I for all sites, i.e. the geomagnetic field do not help sup- 
pressing the electron background for any of the three locations. 
For arrays with energy threshold smaller than ~ 19 GeV (~ 27 
GeV) and O^ax - 40° {0„,ax - 60°), the best site (concerning just 
the suppression of the electron background) is SAC. Again, the 
results obtained for El Leoncito are quite similar to those of 
SAC. For the Namibia site, the suppression of the electron flux 
due to the geomagnetic field begins to be important for energy 
thresholds smaller than ~ 10 GeV, depending on the value of 



Note that the effective area of a given array of lACTs depends 
on the zenith angle of the incident gamma-ray. In order to iso- 
late the dependence of the electron rate with the geomagnetic 
field, this effect is not included in the calculations. Instead, it 
is discussed in section|4]in the context of a particular case, the 
5@5 array. 

As mentioned above, there are regions on the arrival direction 
space for which the positron component is important. Although 
the positron fraction is small at energies close to the £,/, consid- 
ered here, the positron rate can be important in regions of large 
values of £'^'' (small values of E'^"), rising the importance of 
the positron component. FigureQshows contour plots in Aitoff 
projection of Fp„s/Feie for SAC, as a function of azimuth and 
zenith angles, and using the electron flux Jpfl- As expected, 
the positron contribution is more important in the West region 
and for large values of zenith angles. The maximum value of 
Fpos/Feie is 0.82, corresponding to the arrival direction = 60° 
and (p = 270°. In the East region, where takes the small- 
est values and E'^^ the largest ones, FposlFeie is smaller than 
~ 0.1. 

The dependence of F with the energy threshold is also stud- 
ied. Only the arrival directions corresponding to the maximum 
and minimum of F obtained above for £,/, - 3 GeV are used. 
Note that the result for any other arrival direction fall between 



4. Application to the 5@5 array 

The proposed 5@5 array consists of 5 lACTs installed at an 
altitude of 5 km or more [7]. As mentioned above, the site orig- 
inally proposed for the array installation is "Llano de Chaj nan- 
tor", Chile, at an altitude of > 5 km and ~ 200 km northwest 
of SAC. Note that there are several locations at 5 km of alti- 
tude very close to SAC which are suitable for the installation 
of 5@5, near the location of a new astronomical facility under 
consideration 12911. 



The effective detection area of 5@5, for the zenith angle 6 = 
0°, is calculated as 



Aeff{E,e^Q°) = 8.5 



(£/GeV)- 



5.2 



1 + (£/5GeV)* 



ml 



(6) 



The effective area of any array of Cherenkov telescopes changes 
with the zenith angle ifsoll . The distance of the telescope to 
the shower maximum increases, to a good approximation, pro- 
portionally with l/cosfl. As a consequence, the Cherenkov 
photons produced during the shower development illuminate a 
larger area for large zenith angles. Due to this effect the effec- 
tive area increases with increasing by a factor ~ 1 / cos- 6. 
Also, the number of Cherenkov photons that reach the tele- 
scopes decreases with increasing Q as cos^ 9. This causes the 
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Figure 8: F as a function of for tlie three sites under consideration. Only the 
anival directions coiTesponding to the maximum and minimum of F obtained 
for = 3 are used (see text). The electron flux used corresponds to the fit of 
the PAMELA and Fermi LAT data. Solid lines correspond to 6,„„j = 40° and 
dashed lines coiTespond to Omax — 60°. 



energy threshold to increase by a factor ~ 1 / cos^ Q. Therefore, 
the effective area as a function of energy and zenith angle can 
be approximated by llsoll . 



A.//(£,0) 



Aeffi'E COs2 0, = 0°) 

cos^ 6 



(7) 



The left panel of figure |9] shows Agff for the 5@5 array as a 
function of energy for zenith angles = 0°, 30°, 45° and 60°, as 
estimated by combining equations (|6|l and (|7]). 



In order to study the influence of the geomagnetic field on 
the cosmic electron rate measured by 5@5, when it is placed 
in SAC, the parameter F defined by equation Q is also used 
here. The calculation is performed with {F-^ and without (Fq) 
including the zenith angle dependence of the effective area: 



where. 



Ff(e,(p) + F'^'"{9., 



f "" dE T„(E, e, 0) J„{E) A,ff(E,6^0°) 
l^"'" dE J(E)A,ff{E,0^O°) 
dE T,(E, 6, 0) J„{E) A,ff(E, 9) 
l^"'"' dE J(E)Aeff{E,9) 



(8) 
(9) 



(10) 



(11) 



Here a - {ele,pos}. Now the variables in the integrals run 
from zero as the collection area includes the energy threshold. 
Figure [10] shows contour plots in Aitoff" projection of Fq (top 
panel) and F^ (bottom panel), as a function of arriving direction 
of electrons, obtained for the 5@5 array placed in SAC and for 
the electron flux obtained by fitting PAMELA and Fermi LAT 
data (Jpfl)- The values of Fq are larger than those obtained for 
the case studied in the previous section, when a step-function 
effective area with energy threshold £,/, = 3 GeV was assumed 
(table |2]l. Indeed, the maximum and minimum values obtained 
in this case for 5@5 are Fo,,„,„ =0.11 and /^o.moA - 0.31, for 
6»„„„. - 40°. This is due to the fact that for the case of 5@5 the 
product I{E, 6 - 0°) = JeieiE) A^fjiE, 6 - 0°) has a maximum 
at £■ = 5 GeV, as it can be seen from the right panel of figure 
|9] whereas for the ideal case the maximum is reached at £ = 3 
GeV. 

When the dependence of the effective area with zenith angle 
is included, F, changes even more with respect to the ideal case, 
as it is seen from the bottom panel of figure [10] In fact, in the 
region around 6 = 60° and (p - 90°, F~ takes values close to 
one. This can be understood from the behavior of I(E, 0) — 
J(E) AeffiE, 0), shown in figure [9] and E'^'(0, cp) for SAC (see 
top panel of figure O. For (p = 90°, £'^'' decreases from ~ 13 
GeV at 61 = 0° to ~ 8 GeV at 6i = 60°, then the numerator 
of Fr (equation (O) is an increasing function of 0. Also from 
figure [9] it can be seen that the numerator of F. is a decreasing 
function of 0. As a consequence, F~ increases with reaching 
values close to one near = 60°. For <p - 270°, E'j^'^ increases 
from ~ 13 GeV to -25 GeV as goes from 0° to 60°, then, the 
numerator of F~ decreases with 0. However, the denominator 
in equation (O decreases fast enough to make F^. an increasing 
function of the zenith angle. 

Similar results are obtained by using the electron flux 
Jhfl{E). On the other hand, Fq and F~ obtained for "Llano 
de Chajnantor" are slightly smaller than the ones obtained for 
SAC (less than 5% considering all directions) because, as men- 
tioned above, the distance between these two locations is very 
small compared with the scale of variation of the geomagnetic 
field. 

Note that for the case of 5 @5, the importance of the positron 
component decreases when the dependence of the effective area 
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Figure 10: Contour plots of Fq (top panel) and F- (bottom panel) as a function 
of azimuth and zenith angles for the 5@5 aiTay placed in SAC, and using the 
electron flux Jpfl- 

with the zenith angle is included. This is due to the increase of 
the energy threshold of the system with zenith angle. In partic- 
ular, the maximum value reached by F''"'/^^-'' is 0.25, corre- 
sponding to the arrival direction 6 = 47° and <p = 210°. 

5. Conclusions 

In this work we have studied the suppression of the low en- 
ergy region of the cosmic electron flux, due to the effect of 
the geomagnetic field, which is relevant to estimate the back- 



ground for future generations of ground-based gamma-ray de- 
tectors at the lowest energy end of the sensitivity curve. We 
have studied this effect regardless its relative significance com- 
pared to other effects that might be relevant for Cherenkov tele- 
scopes. We have estimated the cosmic electron suppression for 
sites that are, or have been, proposed for the installation of new 
Cherenkov telescopes systems in the southern hemisphere. 

We have considered three sites, two in Argentina (El 
Leoncito and San Antonio de los Cobres - SAC) and one in 
Namibia. To perform these studies we have used numerical 
methods based on the backtracking technique and a multipole 
representation of the geomagnetic field expanded up to order 
10. Considering a step function for the collection area, we have 
found that the largest values of the electron energy, below which 
the flux is dramatically suppressed, correspond to the site SAC. 
The values obtained for El Leoncito are, on average, slightly 
smaller than those corresponding to SAC. We have also found 
that for the Namibia site such values of energy are smaller by a 
factor of two, or even more, depending on the arrival direction 
of the cosmic electron. We found that the total cosmic electron 
flux in SAC is suppressed by 90-97% for an energy threshold 
of 3 GeV, and by 0-75% for an energy threshold of 10 GeV, 
depending on the arrival direction. For an energy threshold of 
19 GeV (27 GeV) none of these locations shields the electron 
background for arrival zenith angles up to 40° (60°). 

We have also studied in more detail the case of the proposed 
5@5 array, including its simulated collection area as a function 
of energy. For this case we have considered two sites: Llano 
de Chajnantor in the Atacama desert, Chile, and SAC in the 
Argentinean Puna, where several nearby locations at 5 km of 
altitude can be chosen. We have found that the geomagnetic 
field effect severely reduces the expected electron rate arriving 
at those locations by approximately 70-90% of the total cosmic 
electrons, depending on the direction of arrival. This suppres- 
sion is slightly larger for the Llano de Chajnantor site, by less 
than 5%. 
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Appendix A. Upper energy cutoff for positrons 

Figure IaTTTI shows contour plots in Aitoff projection of ii^"* 
for SAC, as a function of azimuth and zenith angles. The small- 
est values of Ej^"' correspond to the West region, the opposite 
of what happens for E"^'' . 
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Figure A. 1 1 : Contour plots of E''"' as a function of azimuth and zenith angles 
for SAC. 
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